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Abstract  -  In  13  decerebrate,  artificially  ventilated  cats 
preganglionic  sympathetic  outflow  to  the  heart  was  recorded 
with  ECG  and  ventilation  signal.  A  novel  algorithm  was 
implemented  that  extracts  weighted  events  representing  burst 
occurrences  and  their  size.  A  multiple  threshold  strategy  was 
used  to  separate  bursts.  Weighted  events  yielded  count  signals. 
Spectral  analysis  of  the  count  signal  revealed  a  predominance  of 
a  discharge  synchronous  to  heart  beat  at  cardiac  frequency 
(CF=3  Hz),  a  ventilation  rhythm  at  the  high  frequency  (HF=0.3 
Hz)  of  beat-to-beat  variability,  and  in  9  of  13  cats  a  Mayer's 
wave  lower  frequency  (LF).  The  CF  component  was  —50% 
larger  in  power  than  both  HF  and  also  LF  (when  present). 
Spectral  analysis  at  increasing  count  levels  (i.e.,  with  only  the 
events  with  a  weight  >  the  considered  level)  indicated  that  all 
rhythms  were  carried  by  burst  activity  and  its  modulations.  A 
modulation  index  of  HF  over  CF,  MI  [n/c  f,  was  extracted  from 
dynamic  (i.e.,  ventilation  cycle  by  ventilation  cycle)  folded 
histograms  of  counts.  MIhf/cf  was  significantly  higher  in  the 
group  without  LF  (0.65±0.20,  mean+SD,  n=4)  than  in  that  with 
LF  (0.41±0.07,  n=9).  Burst  activity  can  be  a  key  element  in  the 
interactions  between  cardiovascular  variability  rhythms. 
Keywords  -  Sympathetic  outflow,  burst  activity,  beat-to-beat 
variability,  cardiac  rhythm. 

I.  Introduction 

Rhythms  synchronous  to  cardiovascular  variability  and 
autonomic  control  oscillations  have  been  described  at  various 
levels  of  the  central  nervous  system  and  autonomic  nervous 
system;  i.e.,  the  Mayer's  low  frequency  (LF)  waves  and  the 
respiratory  high  frequency  (HF)  ones  [1,  2,  3].  Nonetheless, 
they  are  accompanied  by  faster  rhythms  that  can  not  be  put  in 
relation  with  cardiovascular  variability  by  means  of  classical 
methods  in  time  and  frequency  domain,  simply  because  they 
are  found  at  a  frequency  equal  or  higher  then  the  cardiac 
frequency  (CF)  itself  [1], 

Complex  interactions  between  rhythms  have  been 
observed  both  in  the  cardiovascular  variability  signals  [4]  and 
in  neural  activity  [5]  and  it  is  well  known  that  non  linear 
mechanisms  do  not  pennit  to  separate  frequency  bands  but  on 
the  contrary  are  capable  to  transpose  the  spectral  contents. 

Multifiber  recordings  hardly  pennit  to  analyze  the 
dynamics  of  single  action  potentials  (spikes),  but  do  offer  a 
rich  infonnation  about  synchronization  patterns  of  different 
fibers  and  their  coupling  to  other  processes  such  as  heart 
contraction  and  ventilation. 

Little  is  known  about  the  processing  of  multifiber  activity, 
beyond  the  classical  methods  of  mediating  over  specific  time 
ranges  (e.g.,  one  cardiac  cycle)  [5]  or  low  pass  filtering. 
Therefore,  it  is  developed  a  novel  algorithm  that  recognizes 
burst  activity  and  represents  it  as  a  series  of  events  weighted 
by  the  estimated  number  of  synchronous  spikes.  Subsequent 
spectral  and  synchronization  analyses  are  carried  out  at 
different  weight  levels  considering  only  bursts  representing  a 
number  of  spikes  equal  or  above  the  specified  level. 


II.  M  ETHODOLOGY 

13  cats  were  decerebrated  under  anesthesia  with  an 
intercollicular  section,  vagotomized  and  artificially  ventilated 
by  a  Harvard  pump  at  18  cycle/min  monitoring  the  pCCh  and 
pCh  levels.  Recordings  were  carried  out  after  cerebral  death 
and  after  the  anesthesia  effect  was  exhausted.  Preganglionic 
sympathetic  outflow  to  the  heart  was  recorded  from  few 
isolated  fibers  of  the  3rd  white  ramus  by  means  of  platinum 
electrodes.  After  amplification  and  antialiasing  filtering  ,the 
neural  signal  was  directly  recorded  and  next  sampled  at  3 
kHz  sampling  rate.  Contemporaneously  the  ECG  tracing, 
ventilation  pump  flow  and  arterial  pressure  were  also 
recorded  and  sampled  at  300  Hz.  Recordings  in  control  (i.e., 
decerebrate)  condition  lasted  about  5  min;  other  conditions 
with  activation  maneuvers  and  denervations  were  also 


Fig.  1.  Neurographic  multifiber  signal  (a).  Rectified  signal 
superimposed  to  multiple  thresholds  (b).  Weighted  event  series:  bar 
height  indicates  the  estimated  number  of  synchronous  spikes  (c). 
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recorded  but  are  not  considered  here.  Sampled  signals  were 
stored  on  a  PC  for  further  processing. 

Multifiber  activity  was  analyzed  by  a  specifically 
designed  program  in  order  to  extract  event  series  in  which 
each  event  was  weighted  by  the  subtended  area  of  the 
rectified  signal  taken  as  the  estimated  number  of  firing  fibers. 
As  shown  in  Fig.la,  the  single  spikes  are  hardly  recognizable 
in  the  neural  tracing  due  to  superpositions  and  noise  but 
bursts  of  activity  are  clearly  detectable. 

The  mean  value  is  first  subtracted  and  the  resulting  signal 
is  rectified.  Fig.  lb  shows  how  the  rectified  signal  is 
compared  with  multiple  thresholds  positioned  at  levels 
increased  by  a  geometric  progression.  Multiple  thresholds  are 
used  to  decide  whether  a  multiple  event  can  be  split  into  two 
subsequent  events  of  less  weight.  The  procedure  is  activated 
when  the  signal  is  above  the  first  threshold  and  works  as 
follows:  1)  the  points  in  which  whatever  threshold  is  passed 
(no  matter  whether  downward  or  upward)  are  detected;  2)  the 
area  subtended  by  the  rectified  signal  between  to  subsequent 
threshold  passages  is  computed;  3)  an  event  is  positioned  at 
the  center  of  this  time  interval  if  and  only  if  the  unitary  area 
estimated  for  a  single  spike  is  exceeded;  4)  the  truncated 
value  of  the  ratio  between  the  area  in  the  interval  and  the 
unitary  area  is  the  weight  of  the  event.  Fig.lc  displays  the 
detected  events  by  means  of  bars  indicating  the  event  times 
by  their  position  and  the  event  weights  by  their  heights. 

It  can  be  noted  that  to  detect  an  event  of  unitary  weight  it 
is  not  sufficient  to  exceed  the  first  threshold  but  also  a 
sufficient  area  must  be  subtended  when  the  signal  is  above 
threshold.  In  addition  when  a  steep  rise  or  fall  is  present  and 
two  subsequent  thresholds  are  passed  in  a  short  time,  no 
event  is  detected  due  to  insufficient  subtended  area.  These 
characteristics  render  the  algorithm  robust  in  respect  to  noise 
and  to  burst  shape  features  and  represent  the  main  rationale 
for  a  multiple  threshold  strategy. 

The  weighted  event  series  is  also  named  counted  signal; 
in  fact,  if  a  pulse  (Dirac's  distribution)  multiplied  by  the 
weight  is  positioned  in  correspondence  to  each  event,  a  signal 
is  obtained  which  carries  the  infonnation  about  the 
instantaneous  firing  rate  exactly  as  if  clusters  of  unitary 
events  were  considered.  Grouping  the  single  events  into 
multiple  events  pennits  to  fix  a  threshold  to  the  weight  of 
events  to  be  considered  and  separate  burst  activity  from 
sparse  activity,  hi  the  following  a  "counted  signal  of  level  n 
will  indicate  that  only  events  with  a  weight  >  n  are 
considered  while  smaller  ones  are  discarded.  All  processing 
procedures  were  applied  to  all  possible  levels  from  level  1 
(i.e.,  the  counted  signal  with  all  events)  up  to  the  maximum 
detected  weight  (usually  15). 

The  spectrum  of  counts  was  obtained  by  a  smoothed 
periodogram  algorithm.  For  a  frequency  analysis  up  to  5  Hz, 
35  sec  long  frames  (66%  overlap)  were  windowed  (Harm's 
window),  FFT  transfonned  and  averaged  over  the  entire  5 
min  recording.  The  same  procedure  was  applied  for  the 
various  count  levels.  As  shown  in  Fig.2  the  spectra  of  counts 
display  the  low  frequency  (LF  around  0.1  Hz)  and  high 
frequency  (HF  =  0.33;  i.e.,  the  ventilation  frequency  in  this 
experimental  condition)  that  can  be  detected  also  in  beat  by 
beat  variabilities  but  also  a  main  spectral  peak  at  cardiac 
frequency  (CF)  that  can  be  analyzed  only  by  processing 
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Fig.  2.  Power  Spectral  Densities  |(spikes  sec)7llz  of  all  events  (level 
1,  see  text)  in  a  case  without  (a)  and  with  (b)  LF.  Same  spectra 
computed  at  level  7  (c  and  d,  respectively). 

directly  the  neural  activity  by  this  or  similar  methods.  The 
Plf>  Phf,  Pcf  power  is  evaluated  by  considering  three  bands 
of  0.01-0.13  Hz,  0.13-0.53  Hz,  and  mean  CF±0.2  Hz, 
respectively.  Both  absolute  and  percent  powers  were 
tabulated. 

Synchronization  to  the  ventilation  cycle  on  a  period  of 
1/HF=3.33  sec  was  evaluated  by  means  of  a  folded  histogram 
of  counts  as  shown  in  Fig.3a.  The  histogram  is  obtained 
simply  by  folding  the  counted  signal  in  correspondence  to 
each  maximum  of  the  ventilation  signal  and  summing  up  the 
weights  of  events  present  at  the  various  phases  of  the 
ventilation  cycle  subdivided  into  bins  of  33.3  msec.  This 
histogram  gives  an  average  and  static  view  of  how  the  neural 
discharge  is  influenced  by  the  ventilation  phase,  hi  order  to 
analyze  how  faster  rhythms  (mainly  CF  rhythm)  are 
modulated  by  ventilation  and  if  they  can  be  phase  locked  to 
it,  a  dynamic  folded  histogram  was  also  considered,  in  which 
the  counted  signal  was  displayed  ventilation  period  by 
ventilation  period,  as  in  Fig.  3b.  This  plot  is  much  like  a  raster 
scan  of  counts  triggered  by  ventilation  but  a  smoothing  is 
obtained  by  grouping  the  events  into  the  33.3  msec  long  bins. 

hi  the  dynamic  folded  histogram  the  activity  synchronous 
with  heart  period  is  clearly  seen  and  also  it  can  be  appreciated 
that  the  different  response  to  ventilation  phase  is  mainly 
obtained  by  means  of  a  modulation  of  the  burst  activity 
synchronous  to  the  heart  beat.  Therefore,  a  modulation  index 
of  HF  over  CF,  MIhf/cb  is  obtained  with  the  following  steps: 

1)  the  value  of  the  highest  and  the  lowest  cardiac  burst  (HCB 
and  LCB  respectively)  is  detected  in  each  ventilation  cycle; 

2)  highest  and  lowest  bursts  are  averaged  over  all  ventilation 
cycles  thus  obtaining  HCBAV  and  LCBAV;  3)  finally, 

MIHf/cf=(HCBAv-LCBAv)/LCBav.  (1) 

III.  Results 

The  average  firing  rate  detected  was  223±64  (mean±SD) 
spikes/sec.  The  activity  was  unequally  subdivided  into  events 
of  different  levels:  e.g.,  events  with  a  weight  of  1  were 
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149±46  /sec,  those  with  weight  5  were  3+0.6/sec  and  with 
weigth  7  were  0.5±0.23/sec. 

Spectral  analysis  of  counts  revealed  that  the  CF  was 
always  the  main  rhythm  in  sympathetic  discharge.  A  clear  HF 
peak  was  always  present;  on  the  contrary  various  levels  of  LF 
were  found  and  in  4  out  of  the  13  cats  this  activity  was 
virtually  absent.  Therefore,  it  was  decided  to  separate  the 
group  into  a  subgroup  with  LF  (SGlr  n=9)  and  into  a 
subgroup  without  LF  (SGlf.not,  n=4). 

hi  SGlr  Pli%=20.91±  15.58,  Phi%=1  9.72+5.9, 

Pcf%=26.2±8.49;  while  in  SGlf.not,  Plf%=1.3+0.56, 
Phi%=22.23±13.12,  Po%=33.08tl7.13. 

hiterestingly,  the  spectral  peaks  are  recognizable  also  in  a 
spectral  analysis  performed  at  level  7,  as  shown  in  Fig.2c,d. 
Indeed,  spectral  peaks  were  detected  even  at  higher  levels, 
thus  indicating  that  burst  activity,  rather  than  sparse  activity, 
carries  the  LF,  HF  and  CF  rhythms.  Only  at  levels  above  12 
spectra  of  counts  became  flat  indicating  that  highest  level 
events  occur  randomly. 

The  folded  histogram  always  revealed  the  effect  of 
ventilation  on  neural  discharge,  hi  Fig.2  the  histogram 
relevant  to  all  events  (level  1)  is  shown  in  black,  while  that 
relevant  to  level  3  (i.e.,  all  events  of  lower  weight  are  cut  off) 
is  represented  by  the  white  bars.  It  can  be  appreciated  that 
also  events  of  higher  weight  are  modulated  by  ventilation 
(this  is  true  also  for  level  up  to  12,  but  it  is  not  shown  due  to 
the  low  scale  of  the  relevant  histograms). 

More  interestingly,  with  the  dynamic  folded  histogram  it 
was  possible  to  detect  the  presence  of  intermittent 
synchronization  patterns  between  the  CF  and  the  HF.  In  the 
example  of  Fig.3b,  it  can  be  appreciated  a  stationary  phase 
relationship  in  the  ventilation  cycles  from  1  to  15  and  from 
52  to  60.  hi  these  epochs  exactly  8  cardiac  cycles  fit  into  a 
ventilation  cycle,  hi  the  intennediate  epoch  (ventilation 
cycles  from  16  to  51)  a  regular  phase  shift  is  on  the  contrary 
present  indicating  quasiperiodicity. 

Moreover,  the  ratio  between  the  HF  and  the  CF  spectral 
components  did  not  display  significant  changes  between 
SGlf  and  SGlf.not  (0.84+0.43  and  0.90±0.70  respectively); 
while,  on  the  contrary,  the  MIhf/cf  obtained  from  the  dynamic 
histogram  was  significantly  higher  in  the  latter  group 
(0.4110.07  and  0.6510.20,  respectively).  This  indicates  that, 
in  absence  of  LF  oscillations,  ventilation  has  a  higher 
modulation  on  the  burst  synchronous  to  the  heart  and  also  this 
effect  is  mainly  non-linear  and,  therefore,  is  not  detected  by 
spectral  analysis. 

v.  Conclusion 

The  proposed  algorithm  of  burst  activity  detection  and 
weighing  extends  the  band  of  observed  rhythms  and  also 
pennits  to  evaluate  its  contribution  in  fonning  slower 
variability  rhythms. 

Spectral  analysis  carried  out  over  higher  burst  levels 
indicates  that  LF  and  HF  variability  oscillations  are  conveyed 
by  burst  activity,  which  in  addition  appears  to  be  synchronous 
to  the  heart  beat  at  the  CF.  The  HF  ventilation  rhythm 
appears  as  a  modulation  of  the  CF  bursts  and  the  relevant 
modulation  index  is  blunted  by  LF  oscillations. 

These  preliminary  results  can  introduce  to  a  novel 
methodology  for  the  analysis  of  rhythm  interactions  in 


Fig.  3.  Cumulative  folded  histogram  of  spikes  in  the  ventilation  cycle 
at  level  1  (a,  black)  and  at  level  3  (a,  white).  Dynamic  (ventilation 
cycle  by  ventilation  cycle)  folded  histogram  (b) 


multiunit  neurographic  recordings. 
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